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ABSTRACT 

We used 3.6, 8.0, 70, 160 /im Spitzer Space Telescope data, James Clerk Maxwell Tele- 
scope HARP-B CO J=(3-2) data, National Radio Astronomy Observatory 12 meter telescope 
CO J=(l-0) data, and Very Large Array Hi data to investigate the relations among PAHs, 
cold (~ 20 K) dust, molecular gas, and atomic gas within NGC 2403, an SABcd galaxy at a 
distance of 3.13 Mpc. The dust surface density is mainly a function of the total (atomic and 
molecular) gas surface density and galactocentric radius. The gas-to-dust ratio monotonically 
increases with radius, varying from ~ 100 in the nucleus to ~ 400 at 5.5 kpc. The slope 
of the gas-to-dust ratio is close to that of the oxygen abundance, suggesting that metallicity 
strongly affects the gas-to-dust ratio within this galaxy. The exponential scale length of the 
radial profile for the CO J=(3-2) emission is statistically identical to the scale length for the 
stellar continuum-subtracted 8 /im (PAH 8 /im) emission. However, CO J=(3-2) and PAH 
8 /im surface brightnesses appear uncorrelated when examining sub-kpc sized regions. 

Key words: galaxies: individual: NGC 2403, galaxies: ISM, galaxies: spiral, infrared: galax- 
ies, radio lines: galaxies 



1 INTRODUCTION 

The James Clerk Maxwell Telescope (JCMT) Nearby Galaxies 
Legacy Survey (NGLS) is an 850 /im and CO J=(3-2) survey 
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with the goals of studying molecular gas and dust in a represen- 
tative sample of galaxies within 25 Mpc. The sample contains 
three subsets: a subset of 31 gal axies from the Spitzer Infrared 
Nearby Galaxies Survey (SINGS: IKennicutt et all (20031) sample; 
an HI flux-limited sample of 36 galaxies in the Virgo Cluster; and 
an HI flux-limited sample of 72 field galaxies. At this point in 
time, CO J=(3-2) observations of the SINGS galaxies and Virgo 
Cluster galaxies have been completed. Aside from the 3.6-160 
data from SINGS that resolve structures less than a kiloparsec in 
size in the closest galaxies, HI Nearby Galaxies Survey (THINGS; 
I Walter et alj|2008h are also available for many of these galaxies. 
Together, these data are a powerful dataset for studying the prop- 
erties of the interstellar medium (ISM) on sub-kiloparsec scales 
within nearby galaxies such as variations in gas-to-dust ratios on 
sub-kiloparsec scales. 

Gas-to-dust ratios within nearby galaxies have been a major 
topic of research since the launch of the Infrared Astronomical 
Satellite (IRAS). Debate on this subject began when dust masses 
derived from IRAS 60 and 100 /im data seemed a factor of 10 too 
low to match what was expected either from the depletion of metals 
from the ISM or from the ratio of gas column density to dust extinc- 
tion (e.g. lYoung et al.ll98dfl98alDevereux & Young|l990|) . More 
recent calculation s of global gas and dust masses using Spitzer 
Space Telescope dWerner et alj|2004h and James Clerk Maxwell 
Telescope data for spiral galaxies produce gas-to-du s t mass ratios 
of ~ 150 (e.g iRegan et alj|2004 iBendo et alj|2006l ; braine et al.l 
l2007h . Although the measured ratios are reliant upon CO to H2 
mass conversion factors that do not vary with metallicity or star 
formation activity, the ratios that were measured are close to what 
is expected from the ratio for the Milky Way derived from t he com- 
parison of gas colum n densities to optical extinction (e.g. [Kriigel 
l2003|j whittet 2003) o r from the depletion of metals from the ISM 
(e.g. lwhittej|2003l ; n i2004l; |Krugelll2008h . Therefore, the natural 
expectation would be that dust emission from kiloparsec-sized re- 
gions should be correlated with the sum of atomic and molecu- 
lar gas emission on those spatial scales. However, since the gas- 
to-dust ratio depends on metallicity and since metallicity varies 
with radius within most galaxies (e.g. Vila-Costas & Edmunds! 
Il992l ; IZaritskv etai]| 19941 ; Ivan Zee et all 19981) . the gas-to-dust ra- 
tio should decrease with radius. For the most part, these assump- 
tions are only beginning to be tested, mainly because the combina- 
tion of HI, CO, and far-infrared dust emission data that is resolved 
on kiloparsec scales has not been available for nearby galaxies. 

Related to this topic is the association between polycyclic 
aromatic hydrocarbon (PAH) emission in the mid-infrared and 
other constituents of the cold ISM. Several observational stud- 
ies using data from either the Infrared Space Observatory (ISO; 
lKessleretai]|l996f) or Spitzer have shown that, in nearby spiral 
galaxies, PAH emission is strongly correlated with cold (~ 20 
K) dust emission in the far-infrared, at least on scales >1 kpc 
(e.g.lMattilaet al.lll999l ; lHaas et alj2002l;lBendo et al.l 2006 ] 2008 ; 



IZhu et al.ll2008l) . Additionally. dRegan et alj|2006h found that PAH 
emission is correlated with CO spectral line emission at submil- 
limetre and millimetre wavelengths within nearby spiral galaxies, 
although these results are based on using radial averages. These re- 
sults strongly suggested that PAH emission is a tracer of both dust 
and gas in the cold interstellar medium (ISM). 

Some of the SINGS/THINGS/JCMT NGLS galaxies have al- 
ready been used in comparisons of HI and CO emission with 
star formation, including measures of star formation that incorpo- 
rate 24 ixm emissio n (e.g. IKennicutt et al] |2~007; Bigie l et al]|2008l ; 
IWilson et al.ll2009h . These studies have found that star formation 



is correlated with molecular gas surface density or with a sum of 
molecular and atomic gas surface density; no correlation exists be- 
tween star formation and HI emission by itself. Since 24 /im dust 
emission is included in these comparisons, the results imply that 
dust emission should be correlated with molecular or total gas sur- 
face densities but not with atomic gas surface densities, although a 
direct comparison between dust, atomic gas, and molecular gas is 
needed. 

As a first look into these topics with the combined SINGS, 
THINGS, and JCMT NGLS data, we compare PAH, cold dust, 
CO J=(3-2), and HI on sub-kiloparsec scales for NGC 2403. This 
SABcd spiral galaxy dde Vaucouleurs et alj|l99ll) w as selected be- 
cause it is relatively nearby dFreedman et alJuOOlL the distance to 
the galaxy is 3.13 ± 0.14 Mpc; ), the galaxy and the dust emis- 
sion is very extended (the optical disc is 21.9 arcmin x 12.3 ar- 
cmin; | de Vaucouleurs et alii 199 lh . While the inclination is 62.9° 
dde Blok et al.l2008f) . the galaxy is close enough to face-on that it is 
possible to distinguish spiral structures and individual star-forming 
regions. As is found in similar late-type spiral galaxies, the regions 
with the stronges t star formation are located outside th e centre of 
of the galaxy (e.g lDrissen et al.ir9 99; Bendo et al. 2008), while the 
cold dust emission peaks in the centre (Bendo et al. 2008). Conse- 
quently, emission related to star formation can be easily disentan- 
gled from emission related to dust surface density. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Mid-infrared data 

We use the 3.6 and 8.0 urn da ta taken with the Infrared Array Cam- 
era (IRAC: lFazio et al.ll2004l) on Spitzer as part of SINGS. The ob- 
servations consisted of a series of 5 arcmin x 5 arcmin individual 
frames taken in a mosaic pattern that covers a 25 arcmin x 25 ar- 
cmin region. The data used for this analysis come from astronom- 
ical observation requests (AORs) 5505792 and 5505536, which 
were separated by four days; these separate observations allow 
for the identification and removal of transient phenomena, particu- 
larly asteroids. The full-width half-maxima (FWHM) of the point 
spread functions (PSFs), as st ated in the Spitzer Observer's Manual 
iSpitz.er Science Center! l2007lfl are 1.7 and 2.0 arcsec at 3.6 and 
8.0 /xm, respectively. Details on the observations ca n be found in 
the do cumentation for the SINGS fifth data delivery (SING S Team] 
l2007lPi 

The final images were created from version 14 of the basic 
calibrated data frames produced by the Spitzer Science Center. 
The data were then processed through the SINGS IRAC pipeline, 
which applies distortion and rotation corrections, adjusts the offsets 
among the data frames, removes cosmic rays, subtracts the back- 
ground, and then combines them together using a drizzle technique . 
A d escription of the tec hnique is presented in IRegan et al.l d2006f) 
and (|SINGS Team 2007). T he flux calibrat ion is expected to be ac- 
curate to 3% iSvitzer Science Centerl l2007). 

The final IRAC images include a residual background with a 
gradient. We interpolated a smoothed version of the background 
outside the optical disc of the galaxy to determine the background 
levels within the optical disc, and we then subtracted this back- 
ground from the data. We then masked out regions contaminated 



http://ssc.spitzer.caltech.edu/documents/som/ 
2 Available at http://data.spitzer.caltech.edu/popular/sings 
/20070410_enhanced_vl/Documents/sings_nfth_delivery _v2.pdf 
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by emission from bright foreground stars (identified as unresolved 
sources in the unconvolved IRAC data with 3.6 fim/8 fim sur- 
face brightness ratios >5). To correct for the diffusion of light 
through the IRAC detector substrate, the 3.6 and 8.0 fim data 
are multipli e d by the "infinite" aperture corrections described by 
Reach et al. I d2005h . which should not be confused with the types 
of aperture corrections used for measuring the flux densities of un- 
resolved sources. The correction factors are 0.944 and 0.737 at 3.6 
and 8.0 fim, respectively. We then subtracted the stellar continuum 
from the 8 fim surface brightness maps (in MJy sr _1 ) using 

I V {PAH 8/im) = I v (8fj,m) - 0.232/„ (3.6^m), (1) 

which was derived bv lHelouetaT] l l2004h . We will henceforth refer 
to this as PAH 8 fim emission so as to distinguish it from the 8 fim 
emission that includes stellar continuum emission. 



2.2 Far-infrared data 



We used the Multiband Imaging Photometer for Spitzer Rie ke et all 
(MIPS; |2004|) 70 and 160 /jm data taken of NGC 2403 as part of 
SINGS. These observations are composed of the AORs 5549568 
and 5549824. These AORs are scan map observations of the galaxy 
that were performed at the medium scan speed. Each scan map is 
~ 0.5° wide and 1° long, which is sufficiently wide to cover the 
entire optical disc of the galaxy. The two AORs were executed three 
days apart so as to allow for the identification and removal of tran- 
sient sources, although most detectable sources usually affect only 
the 24 /im band. The FW HM of the PSFs are 18 arcse c at 70 fim 
and 38 arcsec at 160 ftm iSmtzer Science Centej 2007|). The doc- 
umentation for the SINGS fifth data delivery dSINGS Teamll2007h 
contains additional details on the observing strategy. 

The 70 and 160 /im images were created from ra w data frames 
using the MIPS Data Analysis Tools (MIPS DAT; iGordon et all 
l2005h version 3.10 along with additional processing steps. First, 
ramps were fit to the 70 and 160 fim reads to derive slopes. In this 
step, readout jumps and cosmic ray hits were also removed, and 
an electronic nonlinearity correction was applied. Next, the stim 
flash frames (frames of data in which a calibration light source was 
flashed at the detectors) were used as responsivity corrections. Af- 
ter this, the dark current was subtracted from the data, and an il- 
lumination correction was applied. Following this, short term vari- 
ations in the signal (often referred to as drift) were removed, and 
additional periodic variations in the background related to the stim 
flash cycle were subtracted; this also subtracted the background 
from the data. Next, a robust statistical analysis was applied to 
cospatial pixels from different frames in which statistical outliers 
(which could be pixels affected by cosmic rays) were masked out. 
Once this was done, final mosaics were made using pixel sizes of 
4.5 arcsec 2 for the 70 fim data and 9 arcsec 2 for the 160 fim data. 
The residual backgrounds in the data were measured in regions 
outside the optical discs of the galaxies and subtracted, and then 
flux calibration factors (given as 702 ± 35 MJy sr -1 [MI PS in- 
strumental unit] for the 70 fim data bv lGordonet~ai] ( l2007h and 
41 7±5 MJy sr~ 1 [MIPS instrumental unit] -1 for the 160 ^m data 
bv lStansberrv et aL d2007l) ) were applied to the data. An additional 
70 fim nonlinearity correction given as 

/70Mm(true) = 0.581(/70m™ (measured)) 1 13 (2) 

by iDale et alj (2005) was applied to coarse-scale imaging data 
where the surface brightness exceeded 66 MJy sr~ . 

The final 70 fim image is affected by latent images, which 
are dark images that appear after the array has observed a bright 



source. The latent image effects result in dark streaks to the north- 
west and southeast of the galaxy centre at the edge of the region 
detected at 70 fim (at a =7:37:20 5 =+65:40:00 and a =7:36:35 
8 =+65:32:00). Our algorithms tend to exclude the emission from 
these regions, so the latent image effects only have a very minor 
impact on the analysis. 

We performed several analysis steps in using 70 and 160 fim 
flux densities to calculate dust masses. First, we convolved the 
70 fim data with a kernel that matches the PSF of the data to the 
PSF of the 160 fim data (see Sec. l2.5l below). Next, we regridded 
the 160 fim data so that its coordinate system matched the 70 fim 
data. After this step, we only used data for pixels detected at the 
lOcr level in the convolved data for each wave band so that system- 
atic effects such as offsets in the background would not strongly af- 
fect the results. We proceeded to fit blackbodies B v (T) with tem- 
peratures T modified by A -2 emissivity functions to the 70 and 
160 fim data for each pixel meeting this signal-to-noise criterion so 
as to determine the approximate dust temperatures. Following this, 
we calculated dust masses within individual pixels using 



D 2 f v 



(3) 



K v B v (Ty 

which is a variant of an equation derived bv lffildebrandl dl983h - In 
this equation, D is the distance to the galaxy, /„ is the flux density, 
Kl, represents the absorption opacity of the dust, and B V {T) is the 
blackbody function for the best fitting temperature T. We used the 
160 fim flux densities in Equation [3] because the data were taken 
from the Rayleigh- Jeans side of the thermal dust emission and are 
therefore less sensitive to u ncertainties in temp erature. The value 
for k„ at 160 fim is given bv lLi&Drainel l l200lh as 12.0 cm g 
The dust masses were then converted to dust surface densities by 
dividing the pixels by their surface area in pc 2 . As a test of the ro- 
bustness of the results against systematic offsets in the background, 
we increased the background in one wave band by 3a, decreased 
it in the other by 3a, and then measured the total dust mass and 
the exponential scale length of the radial profile for the dust surface 
density (described in Section|5j. We found that the dust mass only 
varied < 10% and that the scale length varied by only la when this 
was done, even though systematic offsets of this magnitude seem 
highly unrealistic. 

This is a simplistic approach to calculating dust masses, but 
it should provide approximate measurements of the dust masses 
within this galaxy. More sophisti cate d dust models, such as those 
presented bv lDale & Heloul ( 120021) and lDraine et aljj2007l) , include 
multiple thermal components of dust heated by radiation fields with 
variable strengths and also account for the different spectral en- 
ergy distributions of grains of different sizes. Unfortunately, those 
models include many free parameters, and so when only two wave 
bands are available for defining the far-infrared dust SED, as is the 
case for our analysis of subregions in NGC 2403, we cannot accu- 
rately apply complex dust SED mod els. Fortunate l y, wh en the dust 
masses predicted by the models of iDraine et alj d2007l) are com- 
pared with the dust masses calculated using a single thermal com- 
ponent fit to > 70 am data, the masses ag ree to within a factor 
of 2-4 dRegan et alj|2004l ; iBendo et alj|2006h . Therefore, the dust 
masses estimated from a single thermal component should still be 
useful for our analysis. 

It is possible that the observations here may have missed the 
presence of very cold dust (dust with temperatures < 15 K) that 
may not contribute strongly to the 70 and 160 ^m bands but that 
may constitute the bulk of the dust mass within this galaxy. If 
present, this very cold dust emission would be detected at 200- 
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1000 /im. Unfortunately, we do not have any useful data with which 
to search for the presence of such emission. Some archival JCMT 
SCUBA data are available, but nothing is detected with a signal- 
to-noise level of > 3 in the data. Moreover, it is unclear whether 
submillimetre emission that exceeds thermal dust emission mod- 
els fit to < 200 fim originates from large masses of cold dust 
or from dust with unexpectedly high submillimetre emissivities. 
In any case, most current studies on this subject suggest that the 
70 and 160 fim bands can be used without submillimetre data to 
produce reasonably accurate estimates of the dust mass. See Ap- 
pendix|A]for an extended discussion. 



velocity channels. Moment maps were created from the original 
data cube using a mask created from findclumps. We converted the 
CO J=(3-2) data from corrected antenna temperature (T^) to the 
main beam temperature scale by dividing by t\mb = 0.60. The 
CO J=(3-2) emission is not detected beyond radii of ~ 3.75 kpc 
(which corresponds to an angular scale of ~ 250 arcsec along the 
major axis). 

For some of the analysis in this paper, the CO J=(3-2) emis- 
sion is converted into molecular gas surface density. The conversion 
factor is dependent on a comparison between the CO J=(3-2) and 
J=(l-0) line emission. We discuss both of these topics in Section|4] 



2.3 CO data 

2.3.1 CO J =(3-2) data 

The CO J=3-2 line observations were obtained at the JCMT as part 
of the JCMT NGLS, using the 16 element heterodyne array H ARP- 
B with the backend spectrometer ACSIS dBuckle et al]|2009h . The 
FWHM of the beam is 14.5 arcsec. We give an overview of the ob- 
servation and data reduction procedures here. For a more detailed 
explanation of the general JCMT NGLS observing and data reduc- 
tion process for raster maps see Warren et al. (2009, submitted). 

Observations for NGC 2403 by the JCMT NGLS took place 
over two runs in November 2007 and January 2008. We mapped the 
galaxy using a basketweave raster scanning method, with half array 
steps between scan rows/columns. The target area was a rectangular 
region with a position angle of 127° rotated from north through east 
and with a size of 1 1.7 arcmin x 6.2 arcmin, which corresponds to 
half of the optical diameter along both the major and minor axes. 
The correlator was configured to have a bandwidth of ~ 1 GHz and 
a resolution of 0.488 MHz (0.43 kms" 1 at the frequency of the 
12 CO J=3-2 transition) centered at 121 km s _1 . For each scan, we 
integrated for 5 s per pointing within the target field, and scans were 
repeated until we reached the target RMS for the survey (19 mK 
(T\) after binning to 20 km s _1 resolution). After rejecting one 
half of an individual scan that failed the survey quality assurance 
criteria, the total integration time was up to 62.5 seconds per point 
(depending on the number of contributing receptors that were not 
flagged or missing in the array). 

Data reduction and most of the analysis were carried out with 
the version of the Starlink s oftware package that is maintained by 
the Joint Astronomy Centre dCurrieetalj 2008p|. First, data for in- 
dividual receptors with bad baselines were flagged and removed 
from the data. Generally, two to four receptors in a scan were 
flagged. After this, the raw scans were combined into a data cube 
using a sinc(7ra;)sinc(fc7r2;) kernel as the weighting function to de- 
termine the contribution of individual receptors to each pixel. The 
sinc(7ra;)sinc(fc7r2:) function was used instead of other weighting 
functions because it significantly reduces the noise in the maps 
while having a very minor effect on the resolution. The resulting 
cube was then trimmed to remove the leading and trailing scan 
ends outside the target field, where the scan coverage is incomplete. 
A mask was created to identify line-free regions of the data cube, 
and a third-order bas eline was fit to thos e line-free regions. The 
clumpfind algorithm (Willi ams et al.ll994), implemente d as part of 
the Starlink/CUPID task findclumps dBerrv et alj|2007l) . was used 
to identify regions with emission above two times the rms noise 
in a data cube that had been boxcar smoothed by 3 pixels and 25 



Available for download from http://www.jach.hawaii.edu 



2.3.2 CO J=(l-0) data 

We compared the CO J=(3-2) data with CO J=(l-0) data from 
two sources. We primar ily used CO J=(l-0) data taken by 
iThornlev & Wilsorj dl995h using the National Radio Astronomy 
Observatory (NRAO) 12 meter telescope. The observations were 
performed as a series of pointings separated by 30 arcsec that ap- 
proximately cover the central 5.5 arcmin x4 arcmin of the galaxy. 
These data have a FWHM of 54 arcsec and a spectral resolution 
of 2.6 km s _1 . Since the spacing between the pointings is smaller 
than the b eam, we can treat these d ata as a fully sampled map. 
Table 1 by IThornlev & Wilsorj J 19951) gives the integrated line in- 
tensities (based on main beam temperatures) and uncertainties that 
they measured in individual pointings across the optical disc of the 
galaxy; we used these data to construct a CO J=(l-0) image. Up- 
per lim its were treated as for our analysis. See lThornlev & Wilsonl 
dl995h for additional details on the observations and data reduction. 

We also had access to CO J= (l-0) data taken in the BIMA 
Survey of Nearby Galaxies (SONG: iHelfer et al.ll2003h . While the 
6.3 x 5.8 arcsec FWHM of the beam for the data is superior to that 
for both HARP-B and the NRAO 12 m data, the BIMA SONG data 
have problems that make it less suitable for our analysis. A com- 
parison of the BIMA SONG data to the N RAO 12 m data rev ealed 
that only ~ 60% of the flux was recovered ( Heifer et alj2003l) . and 
the area in which sources are detected in the BIMA SONG data 
is significantly smaller than that for either the HARP-B or NRAO 
12 m data. While we still use the BIMA SONG data for a quali- 
tative comparison to other data, we will not use it for quantitative 
comparisons. 



2.4 HI data 

The HI data use d for this analysis were taken by THINGS 
( Walte Tet alj I2008T) . The observations were performed with the 
Very Large Array in the B, C, and D configurations, which ensured 
fairly uniform sensitivity from the largest scales all the way to the 
resolution limit set by the long baselines of B-array. Calibrated, 
continuum-subtracted data cubes were created using the Astronom- 
ical Image Processing System. The uncertainties are dominated by 
the calibration uncertainty, which is 5%. 

THINGS released two data products that differ based on the 
weighting of the data. We used the data with the natural weight- 
ing. In these data, all visibility points in the uv plane are weighted 
equally. Large-scale structures, including the diffuse, extended HI 
emission that may be associated with diffuse interstellar dust, 
should be more easily detected in the data with the natural weight- 
ing. The FWHM of the PSF f or the naturally we ighted data is 
8.75 x 6.75 arcsec. According to Walt er et"aH d2008l) . a comparison 
between the VLA data and si ngle-dish data from the literature com- 
piled bv lPaturel et alj d2003b showed a strong agreement for most 
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Table 1. PSFs Used for Creating Convolution Kernels 



Telescope 


Wavelength 




Source 


/ Instrument 




(arcsec) 




Spitzer 1 IRAC 


3.6 (im 


1.7 


Model (STinyTim) 




8.0 /im 


2.0 


Model (STinyTim) 


Spitzer 1 MIPS 


70 /im 


18 


Empirical 




160 fim 


38 


Empirical 


ICMT / HARP-B 


CO J=(3-2) 


14.5 


Gaussian function 


NRAO 12 m 


CO J=(l-0) 


54 


Gaussian function 


VLA 


HI 


8.75 x 6.75 


Gaussian function" 



a The position angle of the major axis is 25.2° from north to east. 

THINGS galaxies including NGC 2403, indicating that significant 
emission is not missing from th e VLA data becau se of gaps in the 
coverage of the uv plane. See IWalter etai] J2008I) for additional 
details on the observations and data reduction as well as the relia- 
bility of the data. To convert the surface brightnesses S v to atomic 
gas surface densities adust on a pixel-by-pixel basis, we used the 
equation 



M pc_- 
8870 \ 



arcsec^ 
FWHM,, 



Jybeam- 1 kms- 1 ,' 1.1 WIIM 
based on equations given bv lWalteretaT] t2008) 



(4) 



2.5 Convolution kernels 

The various images described above have PSFs with varying 
widths, and the shapes of the PSFs vary as well. In particular, the 
Spitzer PSFs have very pronounced Airy rings whereas most other 
PSFs can be approximated as Gaussian functions. To properly com- 
pare the data to each other, we need convolution kernels that not 
only match the FWHMs of the PSFs to each other but also match 
the shapes to each other. 

We created kernels for th is purpose following the instructions 
given bv lGordon et alj J2008h . We first generated images that rep- 
resented the PSFs in the individual wave bands. For the CO J=(3- 
2), CO J=(l-0), and HI data, we used Gaussian functions. For the 

3.6 and 8 /im data, we used PSFs p roduced w ith STinyTirrQ, a 
PSF simulator designed for Spitzer dKrist| [2002). For the 70 and 
160 /im data, we used empirically-determined PSFs created from 
extragalactic obse rvations of point-like extragalactic sources (see 
Young et al. 2009, for additional information). Table[T]gives a sum- 
mary of these PSFs. 

We then used the equation 



K(x,y) =F~ 



W{w)F[PSF 2 {x,y)] 



F[PSFi{x,y)] 



(5) 



to create convolution kernels K(x, y) that match PSFi to PSF2. F is 
a Fourier transform, W(cu) is a radial Hanning truncation function 
given by 



W(uj) 



i[l+cos(^)] 

<x) > Wq 



(6) 



that is used to suppress high-frequency spatial noise in the resulting 
kernels. For each kernel, we adjusted ljq to produce kernels which, 



when visually inspected, contained as much structure related to 
the PSF transformation as possible while excluding high frequency 
noise. However, for converting from one idealized Gaussian profile 
to another, we created convolution kernels using the equation 



K(r) 



(7) 



where o\ is the width of the narrower Gaussian function that is be- 
ing blurred to match the width 02 of the broader Gaussian function. 

The final kernels that we created convert data to match two dif- 
ferent PSFs. For comparing PAH emission to atomic or molecular 
gas, the kernels match the data to the PSFs of the CO J=(3-2) data. 
For comparing gas and dust emission, the kernels match the data 
to the PSFs of the 160 /im data. For the comparison between the 
CO J=(3-2) and J=(l-0) data, the PSFs of the data were matched 
to the PSF of the NRAO data. Note that we used a pre-made kernel 
created by K. D. Gordon for matching the PSF of the 70 /im data 
to that of the 160 /im data! The convolution kernels were applied 
to the CO J=(3-2) data before spectral line extraction. 



2.6 Abundance information 

Since metallicity has been shown to have a possible effect on the 
CO to H2 conversion factor (discussed further in Section [4} and 
since metallicity variations could explain variations in gas-to-dust 
ratios (discussed further in Section [5}, we need abundance gradi- 
ent information for our analyses. Many authors have measured and 
published gradien ts in 12+log(Q/H) within t his ga l axy, including 
iFierro et all dl986t).IVila-Costas & Edmunds! Jl992h,IZaritskv etall 



1994). iGarnett et all < fl997l) and Ivan Zee et alj Jl998h . After seal 



ing the various gradients so that they correspond to a galaxy 
distance of 3.13 Mpc, we find that the 1 2+log(Q/H) gradie nts 
vary from -0.0774 ± 0.001 4 dex kpc" 1 jyan Zee et al.lll998l) to 
-0.098 ± 0.009 dex kpc" 1 jGarnett et alJl 19971) witha mean of 
-0.084 dex kpc -1 and a standard deviation of 0.009 dex kpc -1 . 
The valu e of 12+log(Q/H) at the centr e of the galaxy is als o 
giv en by I Vila-Costas & Edmunds! dl992h . IZaritskv et all i 19941) . 
and Ivan Zee et alj dl998l) : the mean of their measurements is 8.4 
with a standard deviation of 0.13. We use these values for the abun- 
dance gradients in NGC 2403. Co mpared to the abundance gradi- 
ents for other galaxies rep orted by I Vila-Costas & Edmunds! J 19921) 
and IZaritskv et alj Jl994h . both of whom worked with relatively 
large samples, the abundance gradients measured in NGC 2403 are 
rather typical. 

Although the studies cited above report consistent 
12+log(0/H) gradients for NGC 2403, more recent observa- 
tions of abundance gradients in other galaxies have shown that 
the me thods used in the above references may contain systematic 
errors JPilvu gin 200 1 ; iBresolin et al]|2004l : iRosolowskv & Simon! 
20081 e.g .). Moreover, newer methods for calculating oxygen 
abundances show that the gradients may be significantly shallower. 
However, because no abundance gradient data for NGC 2403 has 
been published using these newer techniques, we will use these 
older abundance gradients for now but discuss how shallower 
abundance gradients could affect the results where appropriate. 



4 Available from 

http://ssc.spitzer.caltech.edu/archanaly/contributed/browse.html 



5 The kernel is available at 

http://dirty.as.arizona.edu/~kgordon/mips/conv_psfs/conv_psfs.html 
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3 IMAGES 



Figure Q] shows in their native resolution the 3.6 ^im, 8.0 /im, 
70 /im, 160 /jm, CO J=(3-2), and HI images for the entire op- 
tical disc of the galaxy. In addition, Figure [2] shows the inner 
12 x 9 arcmin for all of these images as well as the CO J=(l- 
0) images from the NRAO 12 m telescope and from the BIMA 
SONG survey. In interpreting these images, we will refer to the 
HII region marked with the cyan square in Figure [2 ] as VS 44. 
This is region 44 in the catalogue produced by dVeron & Sauvavrg 
1 1965b . and it c orresponds to region 128 i n the atlas of HII re- 
gions mapped by Hodge & Kennicutt ( 1983). VS 44 is the brightest 
source of both Ha emission (Dris sen et al I 1999b and 24 /im emis- 
sion jBendo et alj|2008l) and therefore is the site with the strongest 
star formation activity in the galaxy. Since the region is located 
well outside the nucleus, it can be used to distinguish between ef- 
fects related to radius and effects related to star formation activity 
or infrared surface brightness. 

The stellar emission, which is shown by the 3.6 image, 
looks similar to what is expected for a late-type spiral galaxy; the 
bulge is close to non-existent, and the disc, while well-defined, has 
a clumpy structure. The PAH 8 /im image shows the clumpy dust 
structures in the ISM. Some of this structure is also visible in the 
70 and 160 /im images despite the lower resolution of the data. The 
dust stru cture is fairly typi cal compared to most late-type spiral 
galaxies jBendo et al]|2 007). The CO J=(3-2) and J=(l-0) images 
all appear similar to each other. Some of the structures seen in the 
CO emission, such as the spiral arm structure to the west of the 
centre and a few CO-bright regions like VS 44, are similar to those 
traced by PAH and longer-wavelength dust emission. However, the 
CO maps contain a notable hole close to the centre of the opti- 
cal disc, whereas PAH 8 /im and longer-wavelength dust emission 
is still present in this region. Interestingly, t his hole corresponds t o 
the peak in diffuse X-ray emission found bv lFraternali et alJd2002h . 
We discuss this further in Sections [5] and [6] In contrast to all of 
these images, the HI emission is much more extended; the emis- 
sion is detectable well outside the optical disc of the galaxy. A 
number of shell-like structure s are visible within t he HI data. These 
shells were also identified bv lThilker et al.l l |1998l) , who associated 
them with star-forming regions. While bright objects in many of 
the other wave bands, such as VS 44, do not have HI counterparts, 
some of the shells in the HI image, such as the ones at a =7:37:15 
5 =+65:34:30 and a =7:36:50 5 =+65:37:30, do have counter- 
parts in other wave bands. 

Figure [3] and [4] show the dust surface density and temperature 
maps derived as part of the analysis in Section [2~2l The dust mass 
image exhibits structures that appear similar to those seen at 70 and 
160 /im. However, star-forming regions appear notably enhanced 
in the temperature map. The dust temperatures range from 26 K 
in infrared-bright regions (including VS 44) to 17 K in the fainter 
diffuse regions. Note that the < 20K regions located at approxi- 
mately a =7:37:20 8 =+65:40:00 and a =7:36:35 S =+65:32:00 
may have been affected by latent image effects in the 70 /im band. 
The median temperature is 22 K. These temperatures are similar to 
those measured in other nearby galaxies when single thermal com- 
ponents have been fit to far-infrared ISO or Spitzer data that sam- 
ple the peak of the t h ermal dust emission from 60 to 200 fj,m (e.g . 
Popescu et ail l2002t iBendo et alj |2003l; lEngelbracht et all 12004 
Regan et al. 2004l; lBendo et alJl2006l : |Perez-Gonzalez et al.|[ 2006). 
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Figure 3. Image of the dust surface density in NGC 2403 produced after 
the analysis in Section 12.21 The surface densiti es were corrected for the 
inclination of this galaxy (given as 62.9° by Ide Blok et~ai] <2008l) V The 
image size and orientation is the same as the size and orientation of the 
images in Figure [T] The green circle in the upper left corner shows the 
38 arcsec FWHM of the PSF, which is equivalent to the PSF of the 160 [im 
data. The uncertainties in the data are ~ 20%. The black square indicates 
the location of VS 44. 
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Figure 4. Image of the dust temperatures in NGC 2403 produced after the 
analysis in Section 12.21 The image size and orientation is the same as the 
size and orientation of the images in Figure[T] The green circle in the upper 
left corner shows the 38 arcsec FWHM of the PSF, which is equivalent to 
the PSF of the 160 fj,m data. The uncertainties are ^1 K. The black square 
indicates the location of VS 44. Regions with temperatures of < 20 K at 
a =7:37:20 S =+65:40:00 and a =7:36:35 <5 =+65:32:00 may have been 
strongly affected by latent image effects. 



4 THE CO J=(3-2)/J=(l-0) RATIO AND THE 
CONVERSION FROM LINE INTENSITY TO 
MOLECULAR GAS SURFACE DENSITY 

Given the superior resolution of the HARP-B CO J=(3-2) data 
compared to the older NRAO 12 m CO J=(l-0) data and given the 
superior spatial coverage of the HARP-B data to both the NRAO 
12 m and BIMA SONG data, it would be preferable to use the 
HARP-B data to calculate molecular gas surface densities for com- 
parisons to dust emission, particularly PAH 8 emission. How- 
ever, the conversion from CO line emission to molecular gas mass 
(Xco) is generally based on the CO J=(l-0) line. Therefore, the 
CO J=(3-2)/J=(l-0) ratio needs to be determined to convert the 
CO J=(3-2) line to a molecular gas surface density. 
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Figure 1. Images of the entire optical disc of NGC 2403 in multiple wave bands. The 3.6 fim image traces starlight. The PAH 8 fim image traces primarily 
PAH emission. The 70 and 160 (im images trace cold dust emission. The CO J=(3-2) image shows CO emission associated with molecular gas. The HI image 
shows atomic gas emission. All images are 21 arcmin X 18 arcmin with north up and east to the left. Logarithmic colour scales are used for all images except 
the 1 60 (im image, where the logarithmic scale to the second power was used to enhance the structure in the image, and the CO J=(3-2) image, where a linear 
colour scale was used. The colour scales were selected so as to best enhance the structures visible in the optical disc of the galaxy. Green circles in the upper 
left corners of the images show the beam size when th e FWH M is greater than 10 arcsec; smaller beam sizes are difficult to illustrate in this image. The optical 
disc of the galaxy as defined by de Vaucouleurs et al. 1 1991), is shown as a cyan ellipse in the 3.6 fim image. The blue circles in the PAH 8 fim image show 
the positions of foreground stars that were brighter than or as bright as the continuum emission from the galaxy at 8 fim. These stars were masked out in the 
analysis. The magenta box in the CO J=(3-2) image shows the region observed with HARP-B. Pixels with non-detections in the CO data are set to black. 



For this analysis, we compared the HARP -B CO J=(3-2) data 
to the NRAO 12 m CO J=(l-0) data from iThornlev & Wilson! 

for the reasons discussed in Section [2,3.21 We examined 
both median line i ntensities and intensit i es m easured in the pointing 
positions used by IThornlev & Wilsonl (l995). For these compar- 



isons, we use intensities based on the main beam temperatures. For 
comparing median line intensities, we used CO J=(3-2) data where 
the resolution and pixels are matched to those of the CO J=(l-0) 
data, and we only selected pixels that are 3cr detections in both 
images. The median CO J=(l-0) line intensity in these pixels is 
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Figure 2. Images of the inner 12x9 arcmin of NGC 2403 in multiple wave bands. The 3.6 /im, 8.0 /im, 70 fim, 160 /im, CO J=(3-2), and HI images are 
the same as the images in Figure [T] except that the images are cropped and that the colour scales have been adjusted to show greater co ntrast in the centre . 
The colours in the two additional CO J=(l-0) images (from the BIMA SONG survey and from the NRAO 12 m observations by Thornley & Wilson ( 1995) 
are scaled linearly. The magenta box in the CO J=(l-0) image from the NRAO 12 m shows the region that was covered in the observ ations. Pixels with 
non-d etections in the CO data are set to black. The cyan square ind icates the location of VS 44, which is region 44 in the catalogue by JVeron & Sauvavrel 
[19651) and region 128 in the catalogue by Hodge & Kennicutt 1 1983). See the caption for Figure[T]for additional information on the other symbols and lines in 
the figure. 
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Figure 5. Image of the ratio of the CO J=(3-2)/J=(l-0) line ratio for the 
central 12 X 9 arcmin of NGC 2403. The orientation, size, and central co- 
ordinates of the image have been selected to match those of Figure [2] The 
30 arcsec pixels correspond to the po sitions of the pointed observations per- 
formed bv lThornlev & Wilsotjfl993) . To make this plot, the PSF of the data 
in the CO J=(3-2) data cube were matched to the PSF of the CO J=(l-0) 
data, the CO J=(3-2) line intensities were measured, and then the CO J=(3- 
2) data were rebinned to match the pixel scale of the CO J=(l-0) data. 
The green circle in the upper left comer shows the 54 arcsec PSF of the 
CO J=(l-0) data. The uncertainties in this ratio are mainly dependent on 
the uncertainties in the CO J=(l-0) data; the typical signal-to-noise values 
of the ratios range from 3<j to 8<r. Locations where line emission was not 
measured above the 3cr level and locations affected by artifacts along the 
edge of the region observed in the CO J=(3-2) line were left blank. 

1.56 K km s _1 , while the median CO J=(3-2) intensity is 1.08 K 
km s _1 . This gives a median CO J=(3-2)/J=(l-0) line ratio of 0.69. 
The standard deviation of the ratio across the disc is 0.29. Figure[5] 
shows how the line ratio varies across the centre of the galaxy. The 
ratio does not appear to exhibit any clear structure, although the ra- 
tio is low in the southeast side of the region, which corresponds to 
a location where the CO J=(3-2) ratio is also low. The ratios range 
from 0.11 ±0.3 to 1. 2 ±0.3. 

F or comparison, iThornlev & Wilson] j 1994) and I Wilson et al.l 
dl997h measured ratios for indivi dual giant molecular clo uds in 
M33 ranging from 0.4 to 0.8, and lMauersberger et alj (11999) mea- 
sured CO J=(3-2)/J=(l-0) ratios in 28 nearby spiral galaxies rang- 
ing from 0.2 to 0.7. The median CO J=(3-2)/J=(l-0) ratio for 
NGC 2403 appears high but typical compared to these other mea- 
surements. However, the factor of 10 variation in the CO J=(3- 
2)/J=(l-0) ratio across the disc of NGC 2403 is high compared to 
what was measured by these other surveys, and the minimum and 
maximum fall outside the typical ranges of the objects in the other 
surveys. It is quite possible that the variations in the ratio that we 
observed in NGC 2403 are real and that such variations have not 
been observed in other sur veys either because the the regions were 
unresolved or undetected. IWilson et al.1 d2009h using CO J=(3-2) 
data from the JCMT NGLS found evidence for CO J=(3-2)/J=(l- 
0) ratios that vary by a factor of 10 in NGC 4569 as well. It could be 
possible that further analysis of the ICMT NGLS data will reveal 
that such variations are actually typical for nearby galaxies. 

Given the strong variations in the ratio, we did not feel confi- 
dent about applying a simple scaling term to the CO J=(3-2) data 
to scale it to match the CO J=(l-0) data. We therefore examined 
how the CO J=(3-2)/J=(l-0) ratio could be derived from either the 
radius or the CO J=(3-2) surface brightness. For this analysis, we 
used the CO J=(3-2) and J=(l-0) data that were prepared for Fig- 
ure[5] A fit between log(CO J=(3-2)/J=(l-0)) and a linear function 
of both log(CO J=(3-2)/(K km s~ )) and radius (calculated using 



the inclination of 62.9' circ given bv lde Blok et al.l d2008h ) resulted 
in the lowest reduced \ 2 value. The relation is 

log f \ co J =< 3 ~ 2 > ) = (0.70 ± 0.07) log f J t°. J=(3 ~i 2) ) 
°\'co j=M)j v ' ° V Kkms j ^ 

+ (0.088 ± 0.014) - (°-38 ± 0.03). 

Figure [6] shows the relation between the CO J=(3-2)/J=(l-0) in- 
tensity ratio and CO J=(3-2) intensity, the relation between the 
intensity ratio and radius, and the relation given by Equation [8] 
The dependence of the ratio on CO J=(3-2) intensity is very ap- 
parent, but the dependence of the ratio on radius is fairly weak. 
Nonetheless, including a term for radius in the conversion between 
the CO J=(3-2) intensity and the CO J=(3-2)/J=(l-0) ratio does 
provide a better fit (in terms of the reduced \ 2 ) man a function 
that does not include such a term. We also found that the Pearson's 
correlation coefficient increased from 0.54 for the data in the left 
panel of Figure [6] to 0.81 for the data in the right panel, confirm- 
ing that the extra radial distance term reduced the scatter. Based on 
the scatter around the relation given in Equation [8] the conversion 
from the CO J=(3-2) intensity to the CO J=(l-0) intensity should 
be accurate to within 20%. 

Given these results, we can now scale the CO J=(l-0) to H2 
conversion factor for use with the CO J=(3-2) line. Xco for the 
J=(l-0) transition may depend on several environmental factors. In 
particular, i t has been sh o wn that it may be dependent on metal- 
licitv (e.g. lwilsonlll995l: llsraell Il997l. IzTjOd: iBarone et alj|200d : 
llsrael et alj2003l : IStrong et alj2004llsraell2005h . We have decided 
to use the conversion 

log(X co j=(i-o)) = 12.2 - 2.51og[0/H] (9) 

proposed by llsraell |2000), which is based on determining the 
Xco j=(i-o) needed to keep the ratio of total gas surface den- 
sity to far-infrared surface brightness constant in CO clouds in the 
Large and Small Magellanic Clouds. This equation is also found 
to be in agreement with measurements based on gamma-ray emis- 
sion from the Milky Way ^Strong et al„ 2004). We use this equation 
along with the abundance gradients given in Section l2T6l in one set 
of calculations of the molecular gas surface densities. 

However, recent observational results in which the CO ve- 
locity dispersions, radii, and CO luminosities of giant molecular 
clouds are used to infer their mass and hence Xco j=(i -o) suggest 
that the conversion factor does not vary with metallicity jBlitz et all 
120071 : iBolatto et alJ r2008). contradicting the results derived using 
other methods. Therefore, we also use a second set of molecular gas 
surface density data calculated using Xco j=(i -a) = 1-9 x 10 2 ° 
cm " 2 (K km s" 1 )" 1 , which was derived by IStrong & Mattoxl 
dl996h using models of gamma ray scattering that did not include 
radial varia t ions. 

Israel (1997) indicated that radiation field strength could be 
a second factor that affects X co , although so me authors have no t 
found evidence for such a dependency (e.g. IBolatto et al.l 2008). 
The application of the relation found by Israel ( 1997) to our data 
does not appear straightforw ard. The metric for radiation field in- 
tensity used by llsraell ( 19971) uses the ratio of the far-infrared sur- 
face brightnesses measured using IRAS data to the neutral hydro- 
gen density. Because IRAS data does not sample the Rayleigh- 
Jeans side of the dust SED in the same way that Spitzer data 
does, we would need to convert the far-infrared fluxes derived from 
Spitzer data to match that calcul ated from IRA S data, which may 
not be straightforward. Moreover, llsraell J 19971) derived his relation 
using subregions within dwarf galaxies where radial variations in 
abundances were not present and therefore where variations in the 
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Figure 6. The CO J=(3-2)/J=(l-0) intensity ratio plotted as a function of CO J=(3-2) intensity (left panel), radius (center panel), and the right hand side of 
Equation[8](right panel). The data points are for the 30 arcsec square regions shown in Figure[5] The grey lines are the best fit lines to the data. 



gas-to-dust ratio should not be present, whereas NGC 2403 clearly 
exhibits radial variations in abundances and therefore may exhibit 
radial variations in the gas-to-dust ratio t hat would af fect the ra- 
diation field intensity metric. Fortunately. Ilsraell dl997t) found that 
it should still be possible to derive a reasonable approximation for 
Xco based on its relation to oxygen abundances alone. Therefore, 
for simplicity, we will not attempt to include an additional correc- 
tion to Xco for radiation field intensity, although we still discuss 
the implications of this dependency when appropriate. 



5 COMPARISON OF ATOMIC GAS, MOLECULAR GAS, 
AND DUST SURFACE DENSITIES 

To compare the molecular, atomic, and dust surface densities, we 
first measured the radial profiles of the surface densities as well 
as the gas-to-dust ratios. First, all data were convolved with ker- 
nels that match the PSFs to that of the 160 /im data before radial 
profiles were extracted. The gas-to-dust ratio was calculated on a 
pixel-by-pixel basis from images regridded to match the astrome- 
try of the dust surface density images. We then measured means, 
standard deviations, and uncertainties in mean values within ellip- 
tical annuli with widths along the major axes that were equal to 2 
pixels (or 9 arcsec in the dust and gas-to-dust ratio data, 14.6 arcsec 
in the molecular gas data, and 2 arcsec in the atomic gas data). This 
method for setting the annuli widths was used to avoid problems 
with noise caused by having too few pixels fall within given an- 
nuli, but the resulting annuli are small enough that the data, which 
have a PSF with a FWHM of 38 arcsec, will be Nyquist-sampled. 
Pixels where CO or HI line emission was not detected were set to 
in those data. The ratio of the axes in the elliptical an nuli were set 
to acc ount for the inclination of the galaxy (62.9°; Ide Blok et al.l 
l2008h . 

Figure[7]shows the radial profiles of the HI, H2, and dust sur- 
face densities. Two versions of the H2 surface densities are shown. 
The values shown with a solid line were calculated using a constant 
Xco, and the values shown with the dotted line were calculated 
using the formula for Xco given in Equation [9] For simplicity, we 



will refer to the former as the constant Xco molecular gas surface 
density and the latter as the variable Xco molecular gas surface 
density. Although the standard deviation in the CO J=(3-2) radial 
profiles approaches the same level as the mean values at large radii, 
the uncertainties in the mean values are no wider than the thick- 
nesses of these lines. Hence, the mean values are very well defined 
even though significant variance may be seen in CO intensities at 
specific radii. 

The figure shows that most of the mass in the interstellar 
medium of this galaxy is traced by the HI emission, a result also 
obtained by iThornlev & Wilson! dl995h . The HI surface density is 
relatively flat; it increases slightly towards 4 kpc and then declines 
gradually beyond that radius. The constant Xco molecular gas sur- 
face density decreases exponentially with radius. The scale length 
of the constant Xco molecular gas radial profile is 2.32±0.07 kpc. 
The variable Xco molecular gas surface density, however, is so 
close to flat that it is not possible to accurately fit an exponential 
function to the profile. The dust surface density, like the constant 
Xco molecular gas surface density, decreases with radius. Its ra- 
dial profile has an exponential scale length of 3.16 ± 0.10 kpc. If 
we applied an additional correction for radiation field strength in 
Xco, then this would increase the molecular gas surface density 
in regions with high radiation field intensities, which would corre- 
spond to regions in the centre of NGC 2403 with higher dust tem- 
peratures shown in Figure [4] Therefore, the molecular gas surface 
density radial profile would fall in between those of the constant 
Xco and variable Xco molecular gas surface densities plotted in 
Figure|7] 

Of particular interest is how the radial profiles of the molec- 
ular gas changes when the abundance-dependent correction is in- 
troduced. If the variable Xco molecular gas surface density is an 
accurate depiction of the true surface density, then it implies that 
~ 25% of the ISM at large radii is comprised of molecular gas. 
If true, this would contradict the typical depiction of how the ra- 
tio of molecular to atomic gas varies within spiral galaxies, with 
the centres of galaxies being dominated by molecular gas and the 
outer discs being dominated by atomic gas dWong & Blitzll2002l; 
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Figure 7. Radial profiles of the HI surface density (solid blue line, top 
panel), H2 surface density (green lines, top panel), and dust surface density 
(solid red line, top panel) as well as the total (molecular and atomic) gas- 
to-dust surface density ratio (bottom panel) plotted up to 10 kpc, which is 
equivalent to the edge of the optical disc. For the H2 and gas/dust ratios, the 
solid lines are for quantities calculated using a constant Xco and the dot- 
ted lines are for quantities calculated using a value of Xco that depends on 
12+log(0/H). The standard deviations of the data are plotted as thinner lines 
around all radial profiles except for the H2 surface density and gas-to-dust 
ratios calculated with a variable Xqq, which have standard deviations that 
are similar to the equivalent quantities calculated with the constant Xco ■ 
The uncertainties in the mean values are generally smaller than the thick- 
nesses of these lines. The H2 and dust radial profiles are truncated at the 
maximum radii where the quantities could be reliably measured along the 
major axis (4 and 5.5 kpc, respectively). The gas-to-dust ratio calculated us- 
ing the constant Xco is a l so truncated at the maximum radius at which the 
dust surface densities could be measured. The gas-to-dust ratio calculated 
using the variable Xqq is truncated at 4 kpc to show how the change in 
Xco affects this ratio; beyond this radius, the radial profile merges back 
into the other gas-to-dust radial profile. 



Bigiel et al. 20ol e.g .). However, these prior results rely upon the 
assumption that Xco is constant. This issue should be examined 
further using radial profile data for other galaxies and using deeper 
CO data to determine if molecular gas surface densities calculated 
using an Xco that depends on metallicity are found to still be sig- 
nificant contributors to the total gas surface density at larger radii 
and in other galaxies. If found to be true, then we need to more 
carefully consider how to convert CO intensities into molecular gas 
surface densities and possibly revise the standard depiction of the 
structure of the ISM in spiral galaxies. 

The bottom panel of Figure|7] shows two radial profiles of the 
gas-to-dust ratio based on the two different molecular gas surface 
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Figure 8. Image of the gas-to-dust surface density ratio in NGC 2403. The 
image is 21 arcmin xl8 arcmin with north up and east to the left. Only 
pixels where the dust mass was measured are shown. The molecular gas 
surface densities were calculated assuming that Xco i s constant, although 
the functionality of Xco does not strongly affect the appearance of this 
figure. The green circle in the upper left corner shows the 38 arcsec FWHM 
of the PSF, which is equivalent to the PSF of the 160 fim data. The black 
square indicates the location of VS 44. 



densities used in this analysis. Both radial profiles increase rela- 
tively smoothly with radius. Because the atomic gas is the largest 
component of the ISM, differing assumptions in Xco have only a 
minor impact on the radial profiles within the shown range. Within 
a kpc of the center, the ratio is seen to drop to ~ 100. For compari- 
son, the expected gas-to-dust ratio in the Milky Way near the Earth, 
based either on the measurement of the depletion of metals from the 
gaseous phase of the ISM or on t he ratio of ga s column densi ties 
to op t ical extinctio n, is 100-200 jKrugell [20031 : IWhittedliooH : ILH 
120041 : Kriigel 2008). As we indicated above, the dust masses that 
we have calculated are low compared to more sophisticated models 
of dust emission, so it is unlikely that we have underestimated the 
central gas-to-dust ratio. Instead, it is quite possible that the ISM 
in the centre of NGC 2403 is relatively dust-rich. The gas-to-dust 
ratio reaches a value of ~ 150 at about 2 kpc. At radii of 5.5 kpc, 
which is the limit of where we can accurately measure the radial 
profile of the dust in this galaxy, the ratio reaches ~ 400. 

Figure [8] shows the total gas-to-dust ratio for all pixels in 
which we measured dust masses. To make this image, the data were 
convolved with kernels that match the PSFs to that of the 160 /im 
data (which has a FWHM of 38 arcsec). This image shows that 
the gas-to-dust ratio appears to increase monotonically with radius. 
Most structures visible in Figures Q] and [3] are not at all visible. 
Even VS 44, which has a very high surface brightness in many of 
the images shown in Figure|2] appears indistinct compared to other 
locations at similar radii. The only features that are not radially 
symmetric are the streaks at the northwest and southeast edge of 
the mapped data that are a result of the latent image artefacts in 
the 70 /-im data. The results from Figure [8] suggest that the gas-to- 
dust ratio is primarily dependent on radius and that it does not vary 
significantly between diffuse and star-forming regions. 

Figure[8]is only useful as a qualitative demonstration of how 
the gas-to-dust ratio depends primarily on radius and not on gas 
surface density. As an additional quantitative demonstration of this 
phenomenon, we used techniques similar to what were used by 
lBendoetai1d2008h in their comparisons of PAH 8, 24, and 160 fim 
emission. We used the HI, H2 (calculated using a constant Xco), 
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Figure 9. The dust surface density versus the total (atomic and molecular) 
gas surface density for 45 arcsec 2 regions in NGC 2403. The data were 
extracted from images that had all been convolved with kernels that match 
the PSFs to that of the 160 (im data. Only regions within 5.5 kpc (the limit 
at which the dust masses can be reliably measured) and with S/N ratios 
> 5 were used. Uncertainties in the ordinate are generally smaller than the 
symbols in this plot and the next two plots. See the text in Section |5]for a 
description of how these surface densities were extracted. 



and dust surface density maps that were used to create Figure[8] We 
rebinned these data into 45 arcsec pixels. This pixel size was se- 
lected because it is an integer multiple of the pixel size used for the 
MIPS data that is also close to the 38 arcsec FWHM of the 160 
data. We then extracted the surface densities from each 45 arcsec 2 
region in the data. Because the H I is the dominant component of 
the ISM in this galaxy, the choice of Xco in this analysis has a 
minor impact on the results. 

Figure [9] shows that the total gas density is related to the to- 
tal dust density for these 45 arcsec square regions, but the rela- 
tion exhibits significant scatter. For a given gas surface density, the 
dust surface density could vary by a factor of 5. Figure [Tol shows 
how the gas-to-dust ratio varies with gas surface density. This is an 
alternative form for plotting the data in Figure [9j if a one-to-one 
relation existed between the gas and dust surface densities, then 
the relation in Figure [10] should have a slope of and exhibit lit- 
tle scatter. The data in Figure [10] actually exhibit a negative slope 
(—79 ± 13 (M© pc -2 ) -1 ) and significant scatter. FigurefTTIshows 
how the gas-to-dust ratio varies with galactocentric radius for these 
45 arcsec 2 regions. The relation here shows much less scatter; gas- 
to-dust ratios at a given radius may vary by a factor of 2 or less. 

Using the Spearman's rank correlation coefficient, which is a 
non-parametric correlation coefficient that is useful when dealing 
with data that may exhibit non-linear (i.e. exponential) relations, 
we can demonstrate that the relation in FigurefTTlis more significant 
than the relation in Figure|9] For reference, perfectly correlated data 
would have a coefficient of 1, and uncorrected data would have a 
coefficient of 0. The correlation coefficient for the relation between 
gas and dust surface densities in Figure[9]is 0.73, which would im- 
ply a good correlation. However, the correlation coefficient for the 
relation between the gas-to-dust ratio and radius is 0.91. These re- 
sults demonstrate that the relation between the gas-to-dust ratio and 
radius places stronger constraints on dust surface densities than the 
total gas surface density by itself. Moreover, this quantitative ap- 
proach confirms that the absence of visible structure in Figure[8]is 
not simply a consequence of the poor resolution of the data. Instead, 



Figure 10. The gas-to-dust ratio versus the total (atomic and molecular) 
surface density for 45 arcsec 2 regions in NGC 2403. See the caption for 
Figure|9]for additional information. 
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Figure 11. The gas-to-dust ratio versus radius for 45 arcsec 2 regions in 
NGC 2403. This is can be considered an alternate version of the bottom 
panel in Figure|7] See the caption for Figure|9]for additional information. 



it is a result of the gas-to-dust ratio being primarily dependent on 
radius. However, it is always possible that local variations in the 
gas-to-dust ratio will be visible in higher-resolution data. 

The most obvious driver for this gradient in the gas-to-dust 
ratio is the decrease in metallicity with radius that has been ob- 
served previously (see Section |2~6l >. The gas-to-dust ratio is natu- 
rally expected to increase as abundances decrease simply because 
fewer constituents of dust should be present in the ISM. Gas-to- 
dust ratios significantly higher than the local Milky Way value of 
~ 150 have been measured in many dwa r f galaxies with low m etal- 
licities (e.g. Ilisenfeld & Ferraral Il998l ; IWalter et all 120071) , and 
globally-integrated measurements of gas-to-dust ratios are corre- 
lated with g lobally-integrated 12+log(Q/H) or metallicity measure- 
ments (e.g. Ilssa et alJI 19901 ; Ijames et alj|2002c iDraine et alj|2007l : 
Hirashita et al. 2008). Moreover. Ilisenfeld & Ferraral Jl998h and 



Munoz-Mateos et al.l have found linear relations between 

oxygen abundances and dust-to-gas ratios in nearby galaxies. Al- 
though the dust-to-gas ratios typically decrease faster than the oxy- 
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gen abundances in these relations, the slopes are still on the order 
of 1. 

As stated in 12.61 the typical gradient in 12+log(0/H) that 
has been measured in NGC 2403 is -0.084 ± 0.009 dex kpc" 1 . 
For comparison, we measured the gradient in the logarithm of 
the dust-to-gas ratio (d\og{(Jdust/ o gas ) / dr) to be —0.097 ± 
0.002 dex kpc -1 for NGC 2403. These gradients are so close 
that it implies that the dust-to-gas ratio is primarily dependent on 
metallicity. Alternatively, when we calculated 12+log(0/H) as a 
function of radius, we found the relation between 12+log(0/H) 
and d\og(<7dust / Cgas) /dr had a slope of 0.86 ± 0.02. This slope 
is less than 1, which is c onsi stent with the results obtai ned by 
iLisenfeld & Ferraral i 19981) and lMunoz-Mateos et al] (2009), but it 
is so close to 1 that it implies a linear relationship between the oxy- 
gen abundance and the dust-to-gas ratio. Variations in the ratio of 
oxygen to the other c onstituents of dust could explain the differ- 
ence in the gradients. iGarnett et al. Ul999h . for example, has shown 
that C/O increases systematically with O/H within this galaxy as 
well as within M 101, and since carbon is a primary constituent 
of dust, d\og(cr dust / & gas) / dr would be expected to be steeper 
than the gradient in 12+log(0/H). Additional comparisons between 
radially-averaged dust-to-gas ratios and 12+log(0/H) using all of 
the data from SINGS, THINGS, and the JCMT NGLS would be 
appropriate to determine whether similar results are typically ob- 
tained for all spiral galaxies. 

Some researchers (e.g. Neinin ger et al.l 1 19961 ; llsraell Il997t) 
have suggested that Xco for any galaxy can be determined in a two 
step process. First the gas-to-dust ratio is measured in a region dom- 
inated by atomic gas. Second, that ratio is applied to a region that 
contains CO and dust emission to determine how much molecular 
gas should be present, thus giving Xco ■ However, the results here 
show that the application of this method to spiral galaxies must be 
performed cautiously. The first problem is that this method assumes 
that the gas-to-dust ratio is constant, whereas we have demonstrated 
that it varies with radius. The second problem is that, as explained 
in Section[4] Xco may depend on metallicity, and metallicity will 
vary with radius in most nearby galaxies. Despite these problems, 
it may still be possible to derive Xco within subsections of galax- 
ies by dividing the galaxies into elliptical annuli where both the 
metallicity and the gas-to-dust ratio will remain constant. If regions 
free of molecular gas but containing detectable atomic gas and dust 
emission can be identified within the individual annuli, then those 
locations can be used to measure the gas-to-dust ratio at specific 
radii. Those regions can then be used to infer the gas-to-dust ratio 
for locations at similar radii that contain detectable CO emission 
to determine how much gas should be present and therefore what 
Xco is needed to account for all of the gas. This will then lead 
to a measure of Xco as a function of radius. Unfortunately, it is 
not possible to apply this method to the data used here. To perform 
this analysis correctly, we would need to work with data where the 
PSFs matches the 40 arcsec resolution of the 160 /im data. When 
the CO data are matched to the 160 /im PSF, the area within radii of 
4 kpc (the region in which CO J=(3-2) emission is detected along 
the major axis) contains almost no CO-free regions for comparison 
to regions where CO is detected. 

It is possible that the observations here may have missed the 
presence of very cold dust (dust with temperatures < 15 K) that 
may not contribute strongly to the 70 and 160 bands but that 
may constitute the bulk of the dust mass at larger radii. As argued 
in Appendix lAl it is unclear in cases where submillimetre emission 
has been observed in excess of the ~ 20 K thermal dust models 
fit to the 70 and 160 /im data whether this excess originates from 



large masses of < 15 K dust or very low masses of dust with en- 
hanced submillimetre emissivities. Results from previous analyses 
have suggested that estimates of the dust mass based on the 70 and 
160 /im wave bands sho uld be well within an order of magnitude of 
the expected dust mass dLisenfeld et al.ll2002l; iDumke et al.i r2004 ; 



iRegan et al.ll2004lBendo et alj2006l : lDraine et alj2007l) . Based on 
these previous results, we expect the gas-to-dust ratio shown in Fig- 
ure|7]to be accurate and to not be strongly affected by large masses 
of dust not traced by the 70 and 160 /im bands. It is nonetheless 
possible that the radial increase in our measured gas-to-dust ratio 
could be linked to an increase in the relative fraction of < 15 K dust 
with radius, but to properly determine this will require both submil- 
limetre measurements of the dust emission and models of this dust 
emission that incorporate accurate submillimetre emissivities based 
on laboratory and theoretical results. 



6 COMPARISON OF PAH, CO J=(3-2), AND HI 
EMISSION 

As we are interested in re-examining th e empirical correlati on be- 
tween CO and PAH emission found bv lReganet all ( 2006) using 
the higher-sensitivity HARP-B data, and as they used radial pro- 
files to show this correlation, we will first look at the radial profiles 
for th e CO J=(3-2) inte nsity and PAH 8 /im surface brightness. 
Since IRegan et al.l d2006h used the CO intensities for their analysis 
instead of converting the CO intensities to molecular gas surface 
densities, we will compare the PAH 8 /im surface brightnesses di- 
rectly to the CO J=(3-2) intensities. 

The radial profiles for the HI intensity, CO J=(3-2) intensity, 
and PAH 8 surface brightness are shown in FigureQj] All data 
were convolved with kernels that match the PSFs to that of the CO 
J=(3-2) data (which has a FWHM of 14.5 arcsec) before radial 
profiles were extracted. The techniques used for extracting these 
radial profiles are the same as those used in Section|5] although we 
only used data in which the CO emission was detected to calculate 
the CO J=(3-2)/PAH 8 ratio. The width of the elliptical annuli 
along the major axis in the CO J=(3-2) and CO J=(3-2)/PAH 8 /im 
ratio that we used was still equivalent to 2 pixels (or 14.6 arcsec), 
which does not result in Nyquist-sampled radial profiles. However, 
using narrower elliptical profiles results in excessive noise because 
of the low numbers of pixels falling within the elliptical annuli. 
While the standard deviations in the CO J=(3-2) and PAH 8 /im 
radial profiles are sometimes greater than the plotted radial means, 
just as was the case with the CO radial profiles in Figure [7] the 
standard deviation does not equate to the the uncertainties in the 
means. The uncertainties in the means are actually smaller than 
the thickness of the lines in Figure [T2] except for the ratio of the 
CO J=(3-2) intensity to PAH 8 /im surface brightness. In effect, 
we can actually measure the means to a relatively high precision, 
but we still see significant variance in the data. 

In the top panel of the figure, we can see that the radial pro- 
file of the CO J=(3-2) and PAH 8 /im emission are similar to each 
other. In contrast, the HI radial profile is much broader than either 
the CO J=(3-2) or PAH 8 /im radial profiles. Given the dissimi- 
larity between the PAH 8 /im and HI radial profiles as well as the 
dramatic differences between the appearance of the 8 /im and HI 
images in Figures [T] and [2] we can conclude that the PAH 8 /im 
emission is not associated with HI in this galaxy. We will not dis- 
cuss the HI further but instead focus on the comparison between 
the CO J=(3-2) and PAH 8 /im. 

The CO J=(3-2) and PAH 8 /im radial profiles not only appear 
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Figure 12. Radial profiles of the HI intensity in Jy beam" 1 km s" 1 (blue 
line, top panel), CO J=(3-2) intensity in K km s~ 1 (green line, top panel), 
and PAH 8 fim surface brightness in MJy sr _1 (red line, top panel) as well 
as the ratio of CO J=(3-2) intensity to PAH 8 fim surface brightness (black 
line, bottom panel) plotted up to 10 kpc, which is equivalent to the edge of 
the optical disc. All data were convolved with kernels that match the PSF 
of the CO J=(3-2) data (which has a FWHM of 14.5 arcsec) before radial 
profiles were extracted. The standard deviations of the data are plotted as 
thinner lines around all radial profiles. The uncertainties in the mean val- 
ues are generally smaller than the thicknesses of these lines except in the 
case of the CO J=(3-2)/PAH 8 fim emission ratio, where the uncertainty 
is shown as a shaded grey region. The CO J=(3-2) and PAH 8 fim radial 
profiles are truncated at the radii beyond which the quantities could be re- 
liably measured on the major axis. The ratio in the bottom panel is also 
truncated at the maximum radius where CO J=(3-2) intensities could be 
reliably measured. 



superficially similar but also have similar exponential scale lengths. 
Between 0.5 kpc and 4 kpc, the CO J=(3-2) radial profile has an 
exponential scale length of 2.0 ± 0.2 kpc. The PAH 8 fim emis- 
sion (in the convolved data) has a scale length of 1.99 ± 0.04 kpc 
over the same range. The plot of the ratio of the two in the bottom 
panel of Figure [12] shows that the ratio between 0.5 and 4 kpc has 
a gradient of 0.070 ± 0.017 (Kkrns" 1 ) (MJy sr" 1 ) -1 kpc" 1 , but 
the mean value for the ratio stays between ~ 0.25 and ~ 0.6 at 
radii between 0.5 and 4 kpc. However, while the mean in the CO 
J=(3-2)/PAH 8 fim ratio remains relatively constant over this range 
of radii, the high standard deviation shows that the ratio fluctuates 
greatly at any given radius. The implications of this are discussed 
further below. 



The CO J=(3-2)/PAH 8 fim emission ratio decreases sharply 
in the nucleus of the galaxy. This drop in the ratio corresponds to 
drops in the surface brightness in both wave bands. The drop in 
the ratio may be in part related to variations in the gas-to-dust ratio 
described in Section|5] but that would not explain why the ratio ap- 
pears close to constant at larger radii, nor would it explain why the 
CO and PAH emission also drops near the nucleus. We can exclude 
the possibility that currently-ongoing intense nuclear star formation 
may be inhibiting PAH emission (either through destroying PAHs 
or through changing their ionizati on state) and ionizin g gas in the 
centre of the galaxy ; bot h Ha ( e.g. lDrissen et al.ll 19991) and 24 fim 
images Bend o et all (e.g. r2008h show that the nucleus is relatively 
quiescent compared to other star-forming regions which do have 
PAH 8 fim and CO J=(3-2) counterparts . Because this galaxy is 
classified as having an HII nucleus (e.g. IHo et alJll997T) , and no 
other study indicates the presence of an AGN in NGC 2403, we 
can also exclude the possibility that an AGN could be responsible 
for the decrease in PAH and CO J=(3-2) emission in the nucleus. 

The reason for the decrease in PAH 8 fim and CO J=(3-2) 
in the nucleus of the galaxy may be explained by the simple fact 
that, like many other late-type spiral galaxies, the ISM is clumpy 
and that this galaxy d oes not contain ce ntrally-concentrated gas 
and dust emission (e.g. iBendo et al.|[2007l) . The dumpiness of the 
ISM would also explain why the standard deviations in the PAH 
8 fim and CO J=(3-2) radial profiles in the top panel of Figure [121 
are very large. The CO J=(3-2)/PAH 8 fim emission ratio may de- 
crease towards the centre simply because the two wave bands may 
not be correlated on sub-kpc scales. The comparison of the PAH 
8 fim and CO J=(3-2) images in Section [3] visually demonstrates 
this; some structures in the PAH 8 /im image were not as bright 
as those in the CO J=(3-2) image and vice versa. Additionally, the 
high standard deviation in the radial profile of the CO J=(3-2)/PAH 
8 fim ratio implies that the PAH 8 and CO emission are not 
well-correlated at any given radius. 

For an explicit quantitative analysis to demonstrate that PAH 
8 /im and CO J=(3-2) emission does not trace the same structures 
on small spatial scales, we could not use the same binning pro- 
cedure used to compare the dust surface density versus total gas 
surface density in Section[5] The gaps between detected structures 
in the CO J=(3-2) images made the binned data difficult to inter- 
pret. As an alternative, we used surface brightnesses measured in a 
series of discrete regions. We selected 25 regions with the highest 
intensity in the CO J=(3-2) image and the 25 regions with the high- 
est surface brightnesses in the PAH 8 fim image (convolved with a 
kernel to match its PSF to the PSF of the CO J=(3-2) data). These 
regions were defined using the coordinate system of the CO J=(3- 
2) data; each region is 3 pixels x 3 pixels in the CO J=(3-2) data 
(21.8 arcsec x 21.8 arcsec or 330 pc x 330 pc). Some of the regions 
overlap, and some regions are not centered on peak-like emission 
but instead selected to attempt to cover resolved structures. The 
regions were also selected to avoid locations that may have been 
severely affected by foreground stars. Any locations that were af- 
fected by foreground stars in the PAH 8 fim data were blanked out 
in both data sets before surface brightnesses were extracted. These 
regions are shown in Figure[T3] The extracted surface brightnesses 
are plotted in Figure [141 

The relation between CO J=(3-2) and PAH 8 surface 
brightness differs between the two sets of data, demonstrating the 
presence of regions bright in one band but not the other. The regions 
selected in the PAH 8 fim data exhibit virtually no correlation; the 
Pearson's correlation coefficient for these data is 0.07. The regions 
selected in the CO J=(3-2) data do exhibit a weak correlation; the 
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Figure 13. Images of the PAH 8 fim and CO J=(3-2) emission for the cen- 
tral 12x9 arcmin region in NGC 2403 with green squares overlaid to show 
where intensities or surface brightnesses were measured. The regions se- 
lected by their CO J=(3-2) intensities are shown in the CO J=(3-2) image, 
and the regions selected by their PAH 8 fim surface brightnesses are shown 
in the PAH 8 fim image. See the text for details on the region selection. The 
PAH 8 fim image has been convolved with a kernel to match its PSF to that 
of the CO J=(3-2) image. The green circles in the upper left corner of each 
image show the 14.5 arcsec FWHM of the CO J=(3-2) PSF. Blue regions 
in the PAH 8 fim image show regions strongly affected by stars that were 
excluded from the analysis in both wave bands. 



10.0 



m 
CD 

m 

a> 
a 
as 

Z5 

CO 



CO 

I 
< 

CL 




CO J=(3-2) Intensity (K km/s) 

Figure 14. Comparison of the CO J=(3-2) intensity and PAH 8 fim surface 
brightness for the regions shown in Figure [T3l The regions selected by their 
CO J=(3-2) intensity are shown as green squares. The regions selected by 
their PAH 8 fim surface brightness are shown as red triangles. 



data have a Pearson's correlation coefficient of 0.32. Given that the 
square of the Pearson's correlation coefficient indicates the fraction 
of variation in one variable that can be accounted for by a best fit 
line between the data, we find that, in the CO-selected data, only 
~ 10% of the variation observed in one band can be accounted for 
by a relation with the other band. In the PAH 8/im-selected data, 
< 1% of the variation in one band can be accounted for by a rela- 
tion with the other band. We therefore conclude that, at best, PAH 
8 fim emission is only weakly correlated with CO J=(3-2) emis- 
sion on scales of 330 pc within this galaxy, even though the radial 
profiles have very similar scale lengths. 

We can reject a couple of possible scenarios that would ex- 
plain why the radial profiles for the PAH fim and CO J=(3-2) 
emission should match even while the emission mechanisms are 
not correlated on smaller spatial scales. First of all, PAHs and 
CO do not share a common origin. PAHs are thought to form in 
the atmosp heres of or ou tflows from AGB stars or in interstel- 
lar shocks ((Tielens 2008), while molecular gas is commo nly as- 
sumed to fo rm on the surfaces of interstellar dust grains dKrtigell 
I200il2008l) . It is also unlikely that the radial profiles are matched 
through the Sc hmidt law, wh i ch relates star f ormation to gas 
surface density dSchmidt| [l959; Kenn icutJ[l998l) . While star for- 
mation is strongly dependent on molecular gas surface densities 
in subregions within gal axies jKennicutt et al.ll2007l : iBigiel et al.l 
120081; IWilson et ai]|2009l) . PAH emission is not stro ngly linked to 
other star fo r mation tracers on s u b-kpc scales (e.g. [Calzetti et al.l 
20051 2007: Thi lker et all 120071: iBendo et al.l l2008h . Moreover. 



Bend o et al. I d2008h demonstrated that the ratio of PAH 8 fim emis- 



sion to 24 fim emission (which has been shown to be correlated 
with other star formation tracers) decreases as a function of radius 
within NGC 2403, and the ratio of PAH 8 fim emission to other star 
formation tracers would be expected to vary similarly. 

On possible link between PAH and CO emission is possi- 
ble if molecular cloud formation is triggered in r egion s with stel- 
lar potential wells as suggested by iLerov et al] d2008l) . The new 
molecular clouds would produce CO emission. Additionally, the 
increased starlight in these potential wells (from both the stars orig- 
inally present in the wells and new stars formed from the molecular 
clouds) wou l d heat the diffuse ISM. According to the results from 
Bend o et al. I d2008l) , enhanced PAH emission will be produced by 
locations where the interstellar radiation field has increased. There- 
fore, the stellar potential wells that form molecular gas clouds in 
this scenario would not only enhance CO emission but also en- 
hance PAH emission, thus linking the two on large spatial scales 
and resulting in radial profiles that appear similar for both wave 
bands. However, the enhanced PAH emission does not need to oc- 
cur within the stellar potential wells but could instead occur in re- 
gions surrounding the potential wells. Therefore, regions contain- 
ing molecular gas could exhibit CO emission correlated with PAH 
emission while regions without molecular gas would exhibit PAH 
emission that is uncorrelated with CO emission as shown in Fig- 
ure [14] 

It is also possible that both PAH and CO emission are af- 
fected by similar excitation mechanisms. As stated above, enhance- 
ments in PAH emission may be indicative of increased heating of 
the diffuse ISM, which could be caused by enhancing star forma- 
tion activity (although they still declared that PAH emission was 
not useful for tracing star f ormation on sub-kpc scales). Similarly, 
Knap en & Beckmanl dl996T) postulated based on the comparison of 
CO radial profiles to other star formation tracers that CO could be 
excited by the cosmic rays from star formation, and cosmic rays 
tend to diffuse into the ISM as can be seen by the comparisons of 
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infrared and radio continuum em issi on performed by Murph y" ;t al.l 
d2006^ . lMurphv et al.1 d2006bl) . and iMurphv et al.l J2008I) . There- 
fore, increases in star formation could enhance both PAH and CO 
emission on large spatial scales, which would explain why the ra- 
dial profiles for PAH and CO emission appear to match each other. 

It is also possible that the radial profiles of CO J=(3-2) and 
PAH 8 /im emission are similar simply by coincidence. The radial 
profiles of both normalized by the total gas surface density are ex- 
pected to decrease with radius, although for different reasons. Per- 
haps the multiplicative factors that relate the slope of the CO J=(3- 
2) radial profile to that of total gas surface density just happen to 
match those factors that relate the PAH 8 ^im radial profile to that 
for total g as surface density, at least in this galaxy and the ones 
studied bv lReganetal] d2006h . A more extensive investigation us- 
ing the whole of the JCMT NGLS sample could potentially reveal 
galaxies where the PAH 8 /im and CO J=(3-2) radial profiles are 
not similar (i.e. where the radial profiles do not have similar expo- 
nential scale lengths). If a sufficient number of such galaxies can 
be identified, then it would indicate that what was found here was 
merely a coincidence. 



7 CONCLUSIONS 

We find that the dust surface density depends on the total (atomic 
and molecular) gas surface density and on radius. The gas-to-dust 
ratio varies from ~ 100 in the nucleus to ~ 400 at 5.5 kpc, 
which is the limit at which the dust surface densities can be 
detected. The logarithm of the slope for the gas-to-dust ratio 
(d\og(o dust / a g as) / dr)), which is 0.097 ± 0.002 dex kpc -1 , is 
very close to the slope in 12+log(0/H), which is usually reported as 
between -0.07 to -0.09 dex kpc -1 . This implies that the gas-to-dust 
ratio is strongly dependent on metallicity. Additional abundance 
variations, such as variations in the C/O ratio, could explain the 
small disparity between d\og(o dust /<? gas) / dr and 12+log(0/H). 
Aside from these strong radial variations in the gas-to-dust ratio, 
the ratio does not appear to vary within the galaxy. 

We also find that the PAH 8 fim and CO J=(3-2) radial pro- 
files have statistically identical scale lengths. On spatial scales of 
220 pc, however, PAH 8 /im and CO J=(3-2) emission are un- 
corrected. We propose two mechanisms that could link the PAH 
emission with star formation. First, the CO emission may appear 
associated with PAH emission if molecular cloud formation is trig- 
gered by stellar potential wells. The stars in these potential wells 
would enhance the radiation field that heats the diffuse ISM, which 
would enhance PAH emission. Thus, the PAH emission would ap- 
pear associated with CO emission on large scales but not on small 
scales. Second, the two wave bands may be related through dif- 
ferent excitation mechanisms that are linked through star forma- 
tion and that affect the ISM on large areas. While PAH emission 
does not trace star formation on sub-kpc scales, it is expected to 
be enhanced by increases in the interstellar radiation field that ac- 
company enhanced star formation. Meanwhile, CO emission could 
be enhanced by the increase in cosmic rays that accompanies en- 
hanced star formation. Thus, radially-averaged PAH and CO emis- 
sion would both be linked through star formation, even though PAH 
and CO emission may not be spatially correlated on sub-kpc scales. 
However, we also do not rule out the possibility that the may just 
coincidentally have radial profiles with similar scale lengths. 

This paper is only a first look at the gas-to-dust ratio and the 
relation between PAH and CO emission within the ICMT NGLS 
sample. The SINGS, THINGS, and ICMT NGLS samples were all 



selected to contain many of the same galaxies. The subset of galax- 
ies found in all three samples can be used for an extended compar- 
ison of dust, molecular gas, and atomic gas surface densities and 
for comparisons of PAH to CO emission. Moreover, the gradients 
in the gas-to-dust ratio could then be compared to abundance gra- 
dients to examine the relation between the two further. 
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APPENDIX A: DISCUSSION ON DUST NOT TRACED BY 
THE MIPS DATA 

It is entirely possible that the 70 and 160 /im bands may have 
missed significant masses of very cold dust (dust with tempera- 
tures < 15 K) within this galaxy. However, it is unclear whether 
such dust is present. 

Many surveys working with Spitzer and submillimetre data 
have not found e vidence for very c old d ust. The best examp les are 
the results from bale et al.l d2005l) and iDraine et all d2007l) , who 
worked with Spitzer and submillimetre data for 17 spiral and SO 
galaxies in the SINGS sample and found no submillimetre emis- 
sion in excess of what was predicted by dust models that did not 
include < 15 K dust. Unfortunately, the quality of the submillime- 
tre data used in these studies were qu ite v ariable, so it may not 
have been feasible for lDale et alj ( 120051) and lDraine et ai]( l2007l) to 
detect submillimetre emission in excess of what is expected from 

> 15 K dust within most galaxies. A few groups using very high 
signal-to-noise submillimetre data with ISO or Spitzer data have 
conclusively demonstrated the presence of emission at wavelengths 

> 850 fim in excess of what is expected from ~ 20 K dust withi n 
the Milky W ay Galaxy dReach et al .111 9951; iFinkbeiner et al.lll999f) . 
NGC 4631 dDumke et alj |2004 iBendo et alj 120061). NGC 3310 
dZhu et al.ll2009h. and a few near by dwa rf galaxies i Lisenfeld et al.l 
| 2002l; IGalliano et al l2003t l2005h . While lGalliano et al.ld2003h and 
iGalliano et al.l d2005l) have suggested that this represents emission 
from a high mass of dust at 5-10 K, most other authors have re- 
jected this explanation because it leads to implausibly low gas-to- 
dust ratios. Alternatively, the submillimetre emissivity of dust may 
be enhanced above the A~ 2 emissivity law commonly used in many 
dust models (e.g. iLi & Praind l200ll) . Such dust would include 
grains with exotic shapes, such as fractal dust grains, which would 
have relatively high subm illimetre emissivi t ies relative to their 
absorption cross s ections dReach et al.1 1 19951 ; IPumke et alj [2004 ; 
iBendo et alj|2006t) ; dust in which resonances from i mpurities in 
the dust grains enhances the submillimetre emissivity (Reach et al. 
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ll995l : lBendoetal.ll2006h : or grains with shallow dust emissivities 
(k v oc A"' 3 , with fj < 2) that would only be visible at submil- 
limetre wavelen gths dFinkbeiner et al.lll999HLisenfeld et alj|2002l; 
Zhuetal.|[2009l). Additiona l ly, laboratory studi es l lAgladze et al.l 
19961 : iMermella et all 1 19981: iBoudet et all 120050 and at least one 
theoretical study dMenv et al.l 120070 have shown that dust grains 
exhibit a broad range of emissivities at these wavelength and that 
submillimetre dust emissivities can be strongly enhanced above 
the A -2 emissivity function that is commonly used in astronom- 
ical research. In all of these alternate scenarios, the mass of the 
dust that produces the excess submillimetre emission is relatively 
low, and most of the dust mass should be traced by dust that 
dominates the SEP a t < 450 /xm. Moreover, the resul t s from 



Lisenfeld et ai] feOO20 . iDumke et alj |20040 iRegan et alj ( 120040 . 
Bend o et al.l ( 120060 . and iDraine et al 1 d2007t) have demonstrated 



that, based on the total gas content and expected gas-to-dust ra- 
tios for nearby galaxies, most of the dust mass in nearby galaxies 
can be accounted for by >15 K dust emitting at > 450 /im. 

It is always possible that newer > 160 /im data from either the 
Herschel Space Observatory or SCUBA2 on the JCMT along with 
improved dust models could demonstrate that the use of Spitzer 
data alone will lead to significant underestimates of dust masses. 
However, based on current data and models, the dust masses that 
we measure using only the Spitzer data should provide a reasonably 
accurate estimate of the dust masses. We should not be severely 
underestimating the dust masses because we have not accounted 
for large masses of cold dust emitting > 160 /im. 
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